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Abstract The unexpected development of cycle 24 emphasizes the need for 
a better way to model future solar activity. In this article, we analyze the 
accumulation of spotless days during individual cycles from 1798-2010. The 
analysis shows that spotless days do not disappear abruptly in the transition 
towards an active sun. A comparison with past cycles indicates that the ongoing 
accumulation of spotless days is comparable to that of cycle 5 near the Dalton 
minimum and to that of cycles 12, 14 and 15. It also suggests that the ongoing 
cycle has as much as 20 ± 8 spotless days left, from July 2010, before it reaches 
the next solar maximum. The last spotless day is predicted to be in December 
2012, with an uncertainty of 11 months. This trend may serve as input to the 
solar dynamo theories. 
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1. Introduction 



The sun is a variable star, characterized by the quasi-periodic Schwabe cycle. The 
variation is, among other things, observed by the changing number of sunspots, 
changes in th e total solar irr adiation and by the varying charged particle flux 
from the sun ( Usoskinl . 2008 ). Understanding and predicting this activity is of 
prime importance to the pro tection of sat ellites, astronauts, as well as ground 
based electrical installations ( Jamesl . 120071) . 

The activity is driven by the sola r dynamo, for which we have no complete 
understanding ( Charbonneau . 20051) . Nevertheless, there have been several at- 
tempts to extrapolate future solar activity , based on models and observations 
(jHathawav. Wilson, and Reichmannl . 1 19931) . This also includes the use of spotless 
days as a predictor for fu ture solar activity (|Wilsonl . ll995tlWilson and Hathaway 



2005 . 2006 . 2007 . 20091) . The prediction accuracy varies fro m model to mode l, 



as is clearly demonstrated by the recent extended minimum (|Hathawavll2010l) . 
However, so far no one has used the accumulative number of spotless days to 
predict future solar activity. 

Here we use the Wolf/Zurich/International (WZI) and Hoyt & Schatten (HS) 
sunspot records as a basis for characterizing intervals with low solar activity. The 
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Figure 1. WZI time series (top) and HS time series (bottom). Here blue and green curves show 
the original and interpolated data, respectively. The red curves represent the data smoothed 
with a Gaussian of full width at half maximum (FWHM) of 506 days. The black circles denote 
extrcma points near sunspot minima and maxima. 



cycle minima arc found using a Gaussian filter. Next the accumulation of spotless 
days during individual cycles is analyzed. The analysis includes a condition that 
ensures a stable starting point for the accumulative number of spotless days. 
Based on this condition we compare and discuss the ongoing cycle in relation to 
earlier similar cycles. 



2. The Data 

The sunspot data are available via the National Geophysical Data Center 0. Two 
different tables of sunspot data have been used. The first is a collection of WZI 
sunspot records observed daily since 1818. This collection is quantified by the 
relative sunspot number which is an index that describes the activity of the sun's 
visible disk and is defined as 

R w = K(Wg + f) (1) 

where K is a reduction coefficient that depends on the observational method, 
g is the number of sunspot groups and / is the total number of spots on the 
visible disk. 

The second dataset is a collection of daily sunspot group numbers spanning 
from 1610 to 1995. These data come from 463 independent observers and give a 



more c omplete picture of solar activity on longer timescales (jHovt and Schatten 



1998allbh . From their analysis of historical records Hoyt and Schatten were able 



to construct the sunspot group number and make it compatible with the relative 
sunspot number using a normalization procedure. With this normalization the 



The sunspot data was downloaded on July 15, 2010 at 

http: / /www. ngdc.noaa.gov/stp/solar/ssndata. html 
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index is denned as 



R z = — E k ^ ( 2 ) 

i=i 

where 12.08 is the normalization number which makes the mean i? g 's identical 
with the mean i? w 's for 1874 to 1976 when the Royal Greenwich Observatory 
(RGO) actively made observations of sunspot groups. N is the number of ob- 
servers, fcj is a correction factor, and gi is the number of sunspot groups as 
reported by observer i. The correction factor fcj is defined as 1.000 for the primary 
observer RGO. For all other observers the correction factor is constructed by 
dividing the total number of sunspot groups as reported by the comparison 
observer and by RGO. The details surrounding the creation of the correction 
factor on data befo r e 1874 are somewhat different and is discussed at length by 
Hovt and Schattenl (|l998allbh . 



3. Analysis 



It is important to understand the conditions and limitatio ns of each time series. 
Regarding the WZI time series . lllovt and Schatten (|!998allbT l point out that Wolf 
relied upon correspondents to analyze and send the results to him. The quality 
of these interpretations was sometimes poor because of an unclear definition of a 
group and a sunspot. This is clear from a reexamination of original observations 
from 1819-1833 which shows that some individual sunspots were counted as 
sunspot groups, thus leading to err oneous higher values of the relative sunspot 
number (jHovt and Schattenl . Il995l ). In addition the methodology changed in 
1882 where smaller spots we re counted and the reduction coefficient K changed 
from K = 1 to K = 0.6 (jWaldmeierl . 119611 ). Furthermore, before 1982 the 
index was calculated using a single observer which should give rise to a more 
homogenous time series. As a drawback, all other observations were neglected 
making observational uncertainties difficult to esti mate. Since 19 82 the index 
has been calculated using 20 independent observers ( Usoskin . 20081 ). There is no 
missing observational data after 1849, but in the years between 1818 - 1849 we 
have used linear interpolation to make up for missing data . 

With respect to the HS time series, iHovt and Schatten (|l998al fbh used linear 
interpolation to fill gaps of up to four days for an active Sun and six days for 
a quiet Sun. There is no missing observational data after 1848, but in the years 
between 1798 - 1848 we have used linear interpolation to make up for missing 
data. Moreover, because of large gaps in the time series, the analysis is restricted 
to the period after 1798. 

To locate the time index near sunspot minima and maxima, we use a Gaussian 
filter, with a full width at half maximum (FWHM) of 506 days. This smooths the 
time series enough to locate extrema points near sunspot minima and maxima. 
We then isolate individual cycles based on the located maxima points. The two 
time series, located extrema points and associated cycles are shown in Figure [TJ 
In addition, the intervals of interest, related cycle numbers, associated minima 
and amount of interpolation performed are outlined in Table [1] 
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Table 1. This table illustrates low active intervals based on the zero point condition de- 
scribed in the text. The first column highlight the respective time series, the second is the 
interval of interest, the third display the related cycle number and the forth column show the 
corresponding minimum which was found using the Gaussian filter. TSD is the total number 
of spotless days in the interval and IP is the corresponding interpolation. The last column 
show the estimated accumulated error in the total number of spotless days. 



TS 


Interval 


Cycle 


Minimum 


TSD 


IP (%) 


Error (%) 


WZI 


2007.12.19 - 2010.6.16 


C24 


2009.2.8 


571 










1996.7.4 - 1998.1.9 


C23 


1996.7.21 


152 










1986.6.2 - 1987.7.14 


C22 


1986.6.7 


125 










1976.5.25 - 1977.7.18 


C21 


1976.6.16 


89 










1964.6.5 - 1966.8.10 


C20 


1964.10.8 


168 










1953.10.18 - 1955.10.18 


C19 


1954.6.17 


311 










1943.12.26 - 1945.9.15 


C18 


1944.2.19 


181 










1933.3.13 - 1935.7.29 


C17 


1933.12.11 


390 










1923.1.5 - 1926.7.18 


C16 


1923.6.21 


347 










1912.9.22 - 1916.10.2 


C15 


1913.5.14 


548 










1900.10.31 - 1905.7.28 


C14 


1901.9.20 


639 










1889.8.20 - 1891.12.17 


C13 


1889.10.27 


280 










1877.9.22 - 1883.9.25 


C12 


1878.10.12 


600 










1866.11.4 - 1869.7.14 


Cll 


1867.2.22 


303 










1855.11.3 - 1858.4.5 


C10 


1856.1.27 


380 










1843.6.6 - 1847.7.4 


C9 


1843.9.26 


259 


25.44 


7.797 




1833.12.17 - 1835.6.19 


C8 


1833.12.19 


198 


32.18 


29.86 


HS 


1986.6.2 - 1987.7.14 


C22 


1986.6.5 


125 










1976.5.25 - 1977.7.18 


C21 


1976.7.3 


91 










1964.6.5 - 1966.8.10 


C20 


1964.10.8 


205 










1953.10.18 - 1955.10.18 


C19 


1954.3.3 


349 










1943.12.26 - 1945.9.15 


C18 


1944.2.12 


179 










1933.3.13 - 1935.7.29 


C17 


1933.12.11 


396 










1923.1.5 - 1926.7.18 


C16 


1923.6.21 


339 










1912.9.21 - 1916.10.2 


C15 


1913.5.14 


553 










1900.10.31 - 1903.9.22 


C14 


1901.11.21 


666 










1889.8.21 - 1891.3.27 


C13 


1889.9.23 


269 










1877.9.22 - 1883.5.27 


C12 


1878.10.12 


595 










1866.12.2 - 1868.7.13 


Cll 


1867.2.22 


262 










1855.11.3 - 1858.1.17 


C10 


1856.1.28 


425 










1843.6.6 - 1847.7.30 


C9 


1843.9.23 


270 


2.704 


0.7554 




1833.2.26 - 1835.6.19 


C8 


1833.12.20 


437 


0.8294 


0.4458 




1822.4.12 - 1830.1.24 


C7 


1823.2.5 


982 


1.265 


1.283 




1808.6.10 - 1816.8.27 


C6 


1810.4.10 


1905 


3.732 


0.8849 




1798.4.1 - 1801.6.10 


C5 


1799.1.25 


708 


23.84 


23.85 
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Figure 2. Segment of the WZI time series covering cycle 12 (left) and the accumulated number 
of spotless days during cycle 12 (right). The vertical magenta lines mark the first and last 
spotless day. The corresponding accumulation of spotless days, from the first to the last spotless 
day, is marked by magenta horizontal lines in the right panel. The vertical green line marks the 
zero point spotless day which has AT < 14 days to the previous spotless day. The corresponding 
cutoff is shown by the horizontal green line in the right panel. 



3.1. Determining the Zero Points for Low Activity 



It is of interest to isolate the onset of low solar activity in a given cycle. To 
do this we use the spotless days as a proxy for low activity. The time between 
spotless days varies from months to weeks leading to spotless streaks near solar 
minima. To find the onset of spotless streaks in a given cycle, we isolate all 
the spotless days between the first relevant sunspot maxima and the following 
minima. Next, we calculate the time AT between the neighboring spotless days. 
Finally, we find the last spotless day in this interval, which has AT < 14 days to 
the previous spotless day. This isolated day is then used as the zero point for the 
cumulative summation of spotless days during the relevant cycle. The condition 
AT < 14 days was found such that the zero point lies before the respective 
minima and is close to the onset of the first spotless streak. The procedure is 
illustrated for WZI cycle 12 (C12) in Figure [U where the cutoff in the right 
panel, marks the respective zero point of interest. Repeating the procedure for 
all cycles leads to the final result shown in Figure |H 

3.2. Correspondence between Spotless Days from the HS and WZI Time Series 



It is of interest to compare the ongoing accumulation of spotless days for cycles 
near the Dalton minimum. However, the WZI time series only extends back 
to 1818 and there are gaps of missing data between 1818-1849. Moreover, as 
previousl y mentioned, the rel i abilit y of the WZI time series before 1834 is also in 
question (jHovt and Schattenl . ll995l ). Given that the HS time series ended in 1995 
a direct comparison, to the ongoing accumulation, relies on the correspondence 
between spotless days from the HS and WZI time series. The correspondence is 
best illustrated by plotting the total number of spotless days (TSD), in a given 
cycle, from the WZI time series against similar results from the HS time series. 
The relevant data are outlined in Table [1] and illustrated in Figure [3] This figure 
shows that the scaling of TSD is roughly similar in the two records, except for 
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Figure 3. The total number of spotless days (TSD) during the intervals specified in Table [T] 
from the WZI time series plotted against similar results from the HS time series. Here blue 
show the TSD for cycles C10-C22 with no interpolation and green show the TSD for cycles 
C8-C9 containing interpolation. The scaling is roughly similar in the two records except for 
the outlier of cycle 8. 

cycle 8 which may be due to the uncertain reliability of the WZI time series 
before 1834. Moreover, the change in methodology in the WZI time series, docs 
not seem to affect the number of spotless days after 1882. Thus, the accumulated 
number of spotless days from the two time series are combined in Figure 01 

3.3. Estimating the Accumulated Error 

It is of interest to estimate the accumulated error in the TSD which arise from 
the linear interpolation performed in some intervals. To estimate the error, we 
temporarily remove similar amounts of data in later non-interpolated cycles. 
Next, we employ linear interpolation over this temporary data and calculate 
the TSD using the zero point condition described previously. The result is then 
compared to the non-interpolated TSD in the respective interval. Finally, the 
average deviation is calculated and used to quantify the accumulated error in 
the real interpolated interval. The results are again outlined in Table [1] and are 
also shown in the parenthesis in Figure 2] 
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Figure 4. The accumulative number of spotless days, counted from the first spotless day 
which was found using the zero point condition described in the text. Here red and blue show 
the HS and WZI time series, respectively. The interpolated intervals from both time series is 
shown in green with the estimated accumulated error displayed in the adjacent parenthesis. 
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Figure 5. The remaining number of spotless days, counted from the last spotless day of C24 
in Figure [4] Here red and blue show cycles 12, 14 and 15 from the HS and WZI time series, 
respectively. The estimate on the remaining number of spotless days in C24 is marked in black. 
Note that on the last spotless day of WZI C14 (28th of July 1905) R v = while R g = 1 ± 4. 
Thus, because of this large uncertainty HS C14 is neglected in the estimate of C24. 
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4. Discussion 



The analysis shows that spotless days do not disappear abruptly in the transition 
toward an active sun, but rather there is a gradual decline in the number of 
spotless days. This is also the case for highly active cycles, although the rate 
of decline seems to differ somewhat. Figure 2] shows the ongoing accumulation 
of spotless days compared to previous cycles. It is noteworthy that the ongoing 
accumulation of C24 is comparable to C5 near the Dalton minimum and to cycle 



12,14 and 15 near the modern minimum (jUsoskinl . 120081 ). The cycles following 
C5 and C14 suggest that C25 will be even less active than C24 and this trend 
may serve as input to the solar dynamo theories. 

To estimate the remaining number of spotless days in C24, we shift the zero 
point of Figure S] by 911 days, which corresponds to the present accumulation 
of C24. The result is shown in Figure [5] which illustrates the remaining number 
of spotless days for cycles C12, C14 and C15. There is an obvious difference in 
the last recorded spotless day of C14. From the original data, we note that on 
the last spotless day of WZI C14 (28th of July 1905) R w = while R s = 1 ± 4. 
Thus, because of this large uncertainty, we neglect HS C14 in the following 
calculation. Based on the data related to Figure [S] we find that 20 ± 8 spotless 
days remain, from July 2010, before the next solar maximum. The last spotless 
day is predicted to occur in December 2012 with an uncertainty of 11 months. 
The large uncertainty is related to the difference between C12 and C15. Evidently 
C12 is special in the sense that the last spotless day appears more than 600 days 
later relative to C15. However, C12 is retained in our estimate because of its 
close similarit y to C24 shown in Figure [H 



Previously, IWilson and Hathawavl (120051 ) have shown that the time from the 



first to the last spotless day is strongly correlated (r = 0.949 with a standard 
error of estimate (se) equal to 9.9 months) with the time from first spotless day 
and sunspot minimum occurrence. For C24 the first spotless day occurred in 
January 2004, while a 12-month moving average yields a sunspot minimum in 
December 2008. Thus, the time between the first spotless day and minimum 
occurrence is « 60 months. Therefore, the expected interval of spotless days for 
C24 is 90.7 ± 9.9 months. The time from first spotless day occurrence to July 
2010 is 78 months. Hence, the last spotless day is expected to occur within the 
next 22 months from July 2010, or before May 2012. It sh ould be noted that C12 



i s show n to be an outlier of m 20 months in the analysis of lWilson and Hathaway 



(|2005l ). Thus, there is overlap between the two estimates. 
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